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Resume. 

La methode de calcul de QCD sur reseaux est le seul moyen de calcul non perturbatif 
fonde uniquement sur les bases de QCD. Apres une introduction tres simple aux principes 
de QCD sur reseaux, j'en discute les limites actuelles et la nature des processus qui lui 
sont accessibles. Je presente ensuite quelques resultats frappants dans les domaines des 
■ quarks legers et des quarks lourds. Enfin je tente de deviner les perspectives d'avenir. 

o 
in 
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Abstract 



Lattice QCD is the only non-perturbative method based uniquely on the first principles 
of QCD. After a very simple introduction to the principles of lattice QCD, I discuss its 
Q_i| present limitations and the type of processes it can deal with. Then I present some striking 

results in the light and heavy quarks sectors. Finally I try to guess the prospects. 



X ■ 1 Introduction. 



Let us start with a strong statement. Lattice QCD is within our present knowledge the non- 
perturbative method to solve QCD. It needs no additional assumption beyond QCD, it has 
exactly as many free parameters as QCD itself, that is a coupling constant, or equivalently an 
energy scale (Aq C d for example), and one mass per quark species. 

However, there is a price to pay for this achievement: first, in the proper sense of a price, 
one needs a huge amount of computer time to reach only a limited precision. Second, only the 
simplest processes may be studied by this method, mainly those implying no more than one 
ground state hadron at a time, and in a limited region of the masses and momenta as we shall 
see. Still the improvements are fast, theoretical improvements and numerical ones. The lattice 
community is steadily stepping toward Teraflop machines, dedicated computer are used all over 
the world (but not in France), and we may cherish great hopes. 

The principle of the method is to discretize space-time, and to work in the Lagrangian 
formalism, in a finite volume. Hamiltonian methods with a discretized space and a continuous 
time have also been tried, but they did not achieve by far the same successes and I will skip 
them. 

Several analytic methods are available on the lattice, weak coupling (perturbative) ex- 
pansions, strong coupling expansions, mean field approximations, but numerical simulations a 
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la Monte-Carlo are by far the most fruitful, allowing to fill the gap between the strong and 
weak coupling regimes. 

Numerical simulations on lattice QCD allows you to address almost any question about 
confinement: you may study pure Yang-Mills QCD (i.e. QCD without quarks) and the glueball 
spectrum, you may study at finite temperature the deconfining phase transition as well as the 
chiral one, you may compute the potential between static quarks, the running coupling constant 
starting from non-perturbative data (and this compares quite well with direct perturbative 
estimates), you may learn about the QCD vacuum, its topological properties, its condensates, 
about the theoretical mechanism of confinement (for example you may test the dual Meissner 
effect picture), you may compute the parameters of the chiral Lagrangians, etc. I will here 
concentrate on a few topics directly related to important experimental data: computation of 
masses, structure functions, and electroweak matrix elements. This limited domain is still quite 
wide as we shall see. 

In section ||. I will describe the principles of the lattice calculations, the standard actions, 
the use of path integrals, the Euclidean analytic continuation and the Monte-Carlo method. 
In section |3|. I will summarize and discuss the practical limitations of the method. In section 
f|. I will present some striking and phenomenologically relevant results about light quarks. In 
section ||. I will present the same about heavy quarks. I will then conclude and try to anticipate 
future progress. 



2 The principles and the computational method in short. 
2.1 The pure gauge lattice action. 

Space and time are assumed to be discretized in an hypercubic lattice. It is usual to call a the 
lattice spacing. Let us consider all the oriented links of the lattice. A gauge field configuration 
is described by a set of SU(3) matrices labelled by the oriented links. The relation between 
these SU (3) matrices and the usual continuum gauge field is: 

U^x) = P |e ia9 ° Jo' drAUx+rad)^ j ^ 

where U^x) is the SU(3) matrix attached to the link starting from the site x and going to the 
site x + afi, (/x = 0, 3). go is the bare coupling constant, p, is the unity vector in the positive /i 
direction, i is a color index and Aj, (i = 1, 8) are the so-called "Gell-Mann" Hermitean traceless 
3x3 complex matrices. P means a path ordered product. To the same link, oriented in the 
opposite direction, corresponds the inverse matrix: 

U^(x + afi) = U^(x). (2) 

A gauge transformation is represented by an arbitrary set of 577(3) matrices labelled by 
the sites of the lattice, g(x), and acts as: 

U^x) -> g{x)U i!i {x)g'' 1 {x + aft) (3) 

A plaquette is an elementary square composed by four adjacent links, to which is attached 
the following SU(3) matrix: 

P{x)w = U^(x)U u (x + ap)U- x {x + a^U^ix) (4) 



the trace of which is gauge invariant as it is easy to check. Consequently the simplest pure 
gauge lattice action U is: 



P(*W]} 

(i = l,8). (5) 

for any gauge configuration U, where a gauge configuration is defined as the set U^(x), V(x, n). 

There exists a large number of alternative gauge invariant lattice action which have the 
same limit when a — > as the action in ([|). For a recent summary of the tentative improve- 
ments, see [||]. 



S^] = -]T v,- Re {Tr[l- 
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2.2 Path integrals. 

A field theory is solved once one knows how to compute all Green functions |. The most elegant 
way to compute Green functions, which is used in lattice QCD, is through path integrals. A 
Green function, i.e. the T-product of several operators, is given for example by: 



HdU(x) tlv e iS ^0(x)0'(y) 

< T[0(x)0'(y)) >= = (6) 

/ UdU(x)^e*sM 

X,fJ,,U 

where 0(x) and 0'(y) are operators of the theory, say local polynomials of the elementary 
fields, S[U] is the action defined in (|5]) and dll(x) is a gauge invariant integration measure on 
the SU(3) group. The denominator insures the normalization. In (|6|) there is a T-product of 
two operators, commonly called a two-point correlation function. The generalization to n-point 
correlation functions is straightforward. 

It is important to notice the essential simplification provided by lattice discretization and 
finite volume: the number of integration variables in (||) is finite, while in the continuum theory, 
the path integral is defined in terms of an infinite number of variables with the cardinality of 
the continuum! Still the integral in (^) has a huge number of integration variables: one SU(3) 
matrix per link, i.e. in practice many millions of variables. 

This makes the integral very difficult to compute, all the more so since the exponential 
e iS[{u}} oscillates very fast. The method used to solve the latter problem is to perform an 
Euclidean analytic continuation. 

2.3 Euclidean analytic continuation. 

This consists in continuing the time variable in the complex plane, and using as variables the 
imaginary parts: £4 = ixo,p^ = ipQ. The Euclidean continuation of the two-point correlation 
function in (||) is given by: 



HdUix)^ e- s ^O(x)0'(y) 
' J[dU{x)^ v e-m 

X,fi,U 



(7) 



2 Later on, I will give an example showing how to extract physical quantities from Green function. 



where 

S[U}= £ 4Re{Tr[l-P(x)^]}>0 

X,fl,V 9 

-tiE /^(g;,(x)) 2 >o 

This formula allows a numerical computation as we shall see in the next subsection, thanks 
to the disappearance of the e lS ^ u ^ oscillations. Before turning to that, it is important to stress 
here that once having computed (|7|), an analytic continuation back to the Minkowskian time will 
be necessary. This is straightforward in the simple cases as we shall see, but for systems with 
more than one hadron at the same time the latter analytic continuation far off in the complex 
plane is practically impossible with the presently known methods. This implies that lattice- 
QCD Monte-Carlo simulations are nowadays restricted to physical processes with no more than 
one hadron at one time, except close to the threshold 



2.4 Monte-Carlo. 

The crucial remark is that the positivity of S\U] in (^j) gives to e~ s ^ the meaning of a probabil- 
ity distribution, properly normalized by the denominator in fl7|), and fl7|) is simply the formula 
for a mean value in a probabilistic sense. This suggests an analogy with statistical physics in 
four dimension which has been extensively used to reach a deeper understanding of Euclidean 
QCD. I will skip this and simply notice that although the probabilistic ensemble in (|7|) is the 
huge ensemble of Euclidean gauge configurations, most of them contribute for an exponentially 
suppressed amount. Creutz || suggested to use an algorithm that selects at random gauge con- 
figurations according to the probability law e~ s ^/Z (Z being the denominator in (^)). Such an 
algorithm will discard automatically the configurations whose contribution is negligible. Once 
a large number N of uncorrelated configurations has been produced by the algorithm, the result 
is 

' \{dU{x), v e- s O{x)0'{y) 

HdU(x)^ u e- s N u=i,n WNJ 

X,fl,U 

where O n (x) is the operator 0(x) with the values of the fields in the n th gauge configuration. 
Of course, many different algorithms may be used, improved, tested, and an intense activity is 
going on on this issue. 

Although the computational task is still a formidable one, that needs hundred's of hours 
on supercomputers, it is already a fantastic achievement that this method has proven to be 
practically doable and to lead to non-perturbative predictions that agree quite well with exper- 
iment, within its rather poor numerical accuracy: 0{l/yN). 



2.5 The quark actions. 

Up to now we have considered QCD with only gauge fields. What about quarks ? Quark fields 
are located on lattice sites, and look quite like continuum quark fields except for a different 
normalization: four spin and three color components. The naive action that one is tempted to 
use for quarks on the lattice is: 



+ q(x)mq(x) > (9) 



whose formal limit when a — > is the standard continuum QCD quark action and which is 
gauge invariant as well as chiral invariant when m — 0. 

However, this naive action is not satisfactory because it encounters the "doubling prob- 
lem" . The point is that for one species of quark in (^), the quark spectrum has 16 quarks in the 
continuum limit a — > 0. I have no time to explain why this comes out, let us simply state that 
it has to do with the periodicity of the lattice spectrum. Many proposals have been elaborated 
to solve this problem |J. Wilson's [||] proposal is simple to explain: replace the action (|9]) by 



S qU arks = { ^~ H [<l(. x ) ( r _ 7/J U{x) ll q{x + ap) + (fi -> t7„ -> £//)] + g(x) (m + ~) g(ac) j 

(10) 

The effect of the additional terms proportional to r is to yield a mass 0(r/a) to the 15 
doublers, leaving only one quark with finite mass when a — > 0, as wished. 

Other quark actions that improve the a — > limit are under study ||. It should now be 
clear to the reader that many lattice actions have been proposed or may be proposed for the 
gauge fields and many for the quark fields, leaving aside the many calculational algorithms. 
All these aim at the same theory, QCD. Indeed lattice QCD may be viewed as a class of QCD 
regularization scheme. Discretization cuts off the ultraviolet singularities. The different actions 
are different regularization schemes. Furthermore any regularization procedure of a field theory 
has to be complemented by a renormalization scheme. I will not enter here into a description 
of the subtle question of lattice renormalization. This is also a very ative field of research 0. 

Clearly, all regularization and renormalization schemes must give the same physical quan- 
tities when a — > and the volume goes to infinity. However, due to the numerical uncertainties 
and to the finite value used for a and the volume, these different methods do indeed differ 
somehow in their physical predictions. These differences are taken as a tool to estimate the 
systematic errors due to finite a, finite volume, etc. The statistical errors due to the finite num- 
ber of gauge configurations used in the Monte-Carlo may be estimated from the Monte-Carlo 
itself under some assumptions about the probability distribution of the quantities in question. 
Let us repeat, these different actions and different renormalization schemes are not at all differ- 
ent models, neither do they preclude to the rigor and universality of the calculations performed. 
The are simply different equivalent regularization/renormalization schemes for QCD. 

2.6 An example: computation of a pseudoscalar meson mass and of 

fir' 

Before leaving this section, let us consider in short how one computes a physical quantity on 
a simple example. The axial-current two-point correlation function in Euclidean continuation 
verifies: 



dx < (g(0)7o7 5 g(0)) (q(x, x 4 )7o75g(£, x 4 )) 
= EI <0|g(0) 7 o75<7(0)|n> | 2 

n 

- | <0|g(0) 7 o75g(0)|7T> | 2 



-m n X4 



2m n 

e -m 7T X4 



XA — ' OG 



2m„ 



-/^e— (11) 



where the Euclidean continuation has turned the well known Minkowskian e~ lEx ° time depen- 
dence to the e~ Exi exponential decay. The sum J2 n * s a sum over a complete set in the Hilbert 
space that is coupled to the vacuum by the axial current. At large time, thanks to the expo- 
nential decay, the lightest state dominates the sum. I call the latter state 7r but it can as well 
be a K, a D etc., depending on the mass of the quarks, or even, as we shall extensively use 
later on, a fictitious pseudoscalar meson with quark masses that do not exist in nature. 

The two-point correlation function in ([11]) may be computed by the Monte-Carlo method 
described in subsection |2]4]. Once the calculation performed, the large x 4 exponential slope 
gives the mass m n , and once the latter is measured, the prefactor provides an estimate of the 
leptonic decay constant f n (renormalization introduces here a couple of subtleties which are 
unessential here and that I will skip). An obvious limitation of this method is that the signals 
produced by the excited states are difficult to extract from the statistical noise of the dominant 
ground state: only the ground states are reachable with a suitable statistics, except in the 
NRQCD that will be shortly alluded to in section |5]. 



3 The limitations. 

If it is fair to say that lattice QCD is a rigorous non-perturbative method, it is also necessary 
to insist on its practical limitations nowadays. 

3.1 The quenched approximation. 

The quark action ([TOD is quadratic in the quark fields as in the continuum, and hence, as one can 
read in the text books, the integration of the fermionic variables in the path integral ends up in 
the "fermionic determinant" which is a complicated non-local effective action depending only 
on the gauge fields. In a simpler language, the integration of the fermionic variables amounts 
to compute the dynamical quark loops (i.e. the quark loops generated from the gluons) and 
the result is a non-local interaction between the gluons. On principle there is no obstacle in 
performing this calculation on the lattice and completing to the end the program described in 
the preceeding section. Not only is there no obstacle of principle, but it is currently done by 
several groups. However, these dynamical quark loops increase by several orders of magnitude 
the computation time and increase a lot the difficulty. And this grows even worst for larger 
lattices since the computation time increases like oc n 4 / 3 instead of cx n for a pure gauge theory, 
n being the number of lattice sites. 

It results that for every calculation which demands a rather large number of lattice sites, 
a popular, although unjustified, approximation is performed: the quenched approximation. It 
simply amounts to replace the fermion determinant by 1, or in other words to eliminate all 
dynamical fermion loops such as the following one 




This does not mean that quarks do not propagate. In the example of section 



, the axial 



current creates/annihilates valence quarks. The latter are fully incorporated in the calculation 
which may be pictured as incorporating the whole set of diagrams with any number of gluon 
lines such as this one: 




where the lines at the ends represent the insertion of the axial current. 

The quenched approximation has proven to be quite satisfactory in all processes which 
are dominated by valence quarks. Reversely this approximation has to be avoided whenever 
the sea quarks dominate, and close to the chiral limit when the pion loops (the so-called "chiral 
loops") become important. Indeed it has been shown M that the chiral limit is singular in the 
quenched theory: 



ro, cx m q 
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in the exact theory. 
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(12) 



(13) 



3.2 Limitations on masses and momenta. 

To make it more concrete, the typical values of the lattice spacing practically used are a ~ 
0.05 — 0.1 fm, i.e. a~ l ~ 2 — 4 GeV. The spatial length of the lattice is L ~ 1 — 2 fm; i.e. 
L^ 1 ~ 100 — 200 MeV. To reach meaningful results, the calculation has to be performed with 
a set of parameters such that L^ 1 <C m <C a -1 . Indeed, a -1 is the ultraviolet cut-off of the 
lattice, and L~ l is in some sense an infrared cut-off. A computation performed with masses 
of the order of a -1 contains systematic errors that are totally out of control. Consequently, 
calculations may be performed with quark masses ranging from ~ 50 MeV (jn w ~ 400 MeV) 
to ~ 2.5 GeV. This means that only s and c quarks may be directly considered on a lattice. To 
study quarks u, d, b some kind of extrapolation will be needed. 

The momenta have also to be kept smaller than the ultraviolet cut-off. Furthermore due 
to finite volume and periodic boundary conditions, only a discrete set of values are possible for 
the momenta: 

Pi = — , (2tt/L ~ 0.6 - l.GeV) (14) 

Consequently, for light quarks, no momentum transfer larger than 5 — 10 GeV 2 may be 
studied, and only a few discrete points are measured in the available domain. 



4 The light quarks. 
4.1 The chiral limit. 



We have seen in the preceeding section that it is impossible to put u and d quarks on the 
presently available lattices. One then performs an extrapolation from quarks whose masses 
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Fig. 1: The plot ln(m1/m q ) in the quenched approximation, taken from [|l(| and ||, shows 

the 5 singularity. 



range down as close as possible to the chiral limit. In this process one uses, and simulataneously 
one checks, the chiral theory, i.e. the theory of spontaneous breaking of chiral symmetry. 

For example, in the case of wilson fermions ([K]), one defines massless quarks (the chiral 
limit) as the quarks that combine into a massless pseudoscalar meson, a Goldstone boson 
(the way to measure the pseudoscalar mass has been sketched in section |2.6| ). Notice that 
the parameter m in eq. fllPf ) is the quark bare mass. It is not a zero quark bare mass that 
corresponds to massless renormalized quarks, i.e. to a massless pseudoscalar Goldstone boson. 
Indeed Wilson's regularization procedure fllCf ) breaks chiral invariance and non trivial counter- 
terms are needed to restore chiral symmetry in the renormalization procedure. This strategy, 
of using a Monte-Carlo calculation to fit the bare quark mass that corresponds to the chiral 
limit, is a standard example of non-perturbative renormalization of the mass. 

A recent study of the chiral limit in the quenched approximation |9j exhibits the effect of 
the quenching approximation. Fig. 1 is taken from [10|, ||: it plots as a function of m 



The points to the right of the plot show a nice flat behaviour as expected from chiral theory 
fll3|) , while the points to the left show the effect of the 5 singularity due to quenching (|T^) . One 
learns from this plot that quenching becomes really misleading for quark masses below about 
one half of m s . On the other hand, the linear extrapolation from masses above m s /2 seems 
to be compatible with the chiral theory, and in fair agreement with experiment: paradoxically, 
within the quenched approximation, it seems better to study the chiral limit by extrapolating 
from not too light quarks. This restricts of course the expected accuracy. 

4.2 f n /m p 

One physical scale is necessary to renormalize lattice QCD in the chiral limit. This is necessary 
to fix the scale of the lattice spacing a. Indeed, in massless QCD there is one free scale, usually 
referred to as Aq C d in perturbation theory. Hence, lattice calculations can only predict ratios 
of chiral quantities. For example the ratio j^jvn p (in practice u and d quarks can be considered 
as massless within the present accuracy) has been studied |ll|] in the quenched approximation 
with high statistics, and the results are reported in table 1. 



decay 


finite volume 


inf. volume 


exp. 


U/mp 


0.106(9) 


0.106(14) 


0.121 


fa/nip 


0.121(6) 


0.123(9) 


0.148 


Fp/m p 


0.177(21) 


0.173(29) 


0.199 




0.217(19) 


0.253(35) 


0.219 



Table 1: Light quarks decay constants scaled to the p mass. 



As can be seen, within errors, the quenched lattice agree with experiment (except for fx) 
with a tendency to lie below experiment (except for F^). Whether this tendency is a systematic 
quenching effect is not yet known. 

4.3 mjsr/mp. 

The m^jm p is another ratio concerning light quarks that lattice compute and that can be 
compared with experiment. The first calculations have found a ratio around 1.4-1.5 against an 
experimental number of 1.20. It is now understood that this was simply due to the fact that 
the quark masses where too heavy in the calculations, for the reasons explained in the section 

The High Energy Monte-Carlo Grand Challenge (HEMCG) has performed an unquenched 



calculation |12| , i.e. including the quark loops, oim N /m p as a function of m^/mp, that we report 
in fig. 2. There is an indication that /m p does indeed decrease when m w /mp decreases, and 
the lattice results seem to point toward the experimental point represented on the plot by a 
question mark. The quenched results fl3| are not significantly different from the unquenched 



ones, they confirm the agreement with experiment down to lower values of m n /m p . 
4.4 Structure functions. 

Structure functions are important hadronic quantities about which lattice may say a word. One 
usually defines the n th moment by: 

f dxx n F 1 (x, q 2 ) = -L I dxx n ~ 1 F 2 (x, q 2 ) =< x n > (15) 

the equality being valid for q 2 — ► oo. 

Lattice calculations have been performed for the lowest momenta some time ago [14] and 
seem to be confirmed by recent ones |L5[ which use lower quark masses. The latter find: 



< x u >= 0.42(4); [exp : 0.28] 

< x d >= 0.18(2); [exp: 0.11] 

< x 2 u >= 0.12(2); [exp : 0.08] 
< x 2 d >= 0.050(8); [exp : 0.03] 



where [exp] refers to the valence quark momenta fitted by fll~6| . 
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Fig. 2: m^/m p as a function of m^/rrip, as computed with dynamical quarks (unquenched) 

by the HEMCG collaboration. 



It is not clear to me whether the quenching approximation may be blamed for this dis- 
crepancy. 



5 The heavy quarks. 

We have seen very few of the many lattice results about light quark physics, the preferred field 
for the project ELFE. Let us now turn to the realm of a r-charm factory, namely heavy flavor 
physics. Lattice calculations have provided a lot of results concerning masses, leptonic decay 
constants and semileptonic decay form factors. Notice that nonleptonic decays are out of reach 
because they include more than one hadron in the final state H. The spectroscopy of heavy 



quarkonia is also extensively studied including some excitated states [17] . 



5.1 The charm. 



Table 2, taken from [pl|, reports the leptonic decay constant for charmed mesons. As can 
be seen, the different lattice predictions agree quite well with one another as well as with 
experiment within the large present experimental errors. Notice the unquenched HEMCGC 
results, which also agree with the quenched ones within their expected larger errors. As may 
be seen, the lattice have an accuracy of around 10-15 % for fr> and fo s which were predicted 
at a time where no experimental data was available. In the high statistics calculations, the 
statistical error is as small as 5 %, the systematic errors being larger, as estimated by the 
authors themselves, or as indicated by the discrepancy between different groups^. 

3 Differcnt analysis methods may be used, that should agree in the continuum limit. Their difference for finite 
lattice spacing and finite volume gives an estimate of the systematic errors. This comparison is sometimes done 
inside one collaboration, but it can also be extracted from the comparison of different collaboration's results. 





f D (MeV) 


fo. (MeV) 


/d s / Id 


BLS 


208 ± 9 ± 37 


230 ± 7 ± 35 


1.11 ±.02 ±.05 


UKQCD 


185l|lf 


212l|lf 


1.18 ±0.02 


PWCD 


170 ± 30 




1.09 ±.02 ±.05 


APE 


218 ±9 


240 ±9 


1.11 ± .01 


ELC 


210 ± 15 


227 ± 15 


1.08 ± .02 


LANL 


241 ± 19 


266 ± 15 




HEMCGC 


215±5±53 


288 ± 5 ± 63 




Exp. 








WA75 




232 ± 45 ± 52 




CLEO 




344 ± 37 ± 67 




BES 




434111H 





Table 2: Dt s -j leptonic decay constants from different groups, compared to experiment (in this 
convention: f n = 132 MeV). 





A(o) 


U(0) 


Mo) 


A 2 (0) 


BKS 
UKQCD 
APE 
ELC 
LANL 


0.90(8)(21) 
0.67l;° 8 7 
0.72(9) 

0.60(15)(7) 
0.73(6) 


1.43(45) (49) 
0.981;}° 
1.00(20) 
0.86(24) 
1.24(8) 


0.83(14)(28) 
0.701;? 7 
0.64(11) 
0.64(16) 
0.66(3) 


0.59(14)(24) 
0.681;^ 
0.46(34) 

0.40(28)(4) 
0.45(19) 


Exp. 

(average) 


0.77(4) 


1.16(16) 


0.61(5) 


0.45(9) 



Table 3: D — > semileptonic form factors at q 2 = from different lattice groups, compared 
to experiment. 




Fig. 3: fpMp 2 from pOfl . The two points at the extreme left are static quarks. The points to 

the right correspond to moving quarks. They are in the region of the D meson. All points 
agree. The vertical line corresponds to the B meson. The slope of the fit indicates the size of 

non-leading, 0(1/ Mp), corrections to HQET. 



In table 3, the semileptonic form factors for D decays are quoted again from |Tj|. The 
experimental average are quoted from [|l^]. Again the comparison with experiment is quite 



satisfactory for /+ and A\, satisfactory for V, while A2 does not show a discrepancy within 
larger experimental and theoretical errors. As we have already mentioned, the finite volume 
lattices can only give a few points in momentum space. Practically, for D semileptonic decay, 
in most cases only two points are in the physical q 2 domain. Luckily, however, the masses of 
D and K are such that these two points turn out to be at q^aximum anc ^ c l° se to g 2 = 0, that 
is they bracket the physical domain, thus limiting additional errors that would have resulted 
from an extrapolation. 

It is now clear that lattice methods are in a good position to give their best in this domain 
of D and D s (semi)leptonic decays. Indeed, as already stressed, c and s quarks are in the mass 
range that is easy to deal with on the lattice, and we have just seen that the q 2 physical domain 
is well accessed. It is useful to stress that a better experimental accuracy for these numbers, 
and particularly leptonic decay constants, from a r-charm factory or any other source, would 
be extremely useful to check really the precision achieved by lattice methods. 

5.2 Extrapolation to beauty. 

The b quark is too heavy to be directly studied on present lattices, since m& > a -1 . It is thus 
necessary to resort to an extrapolation. Similarly to the light quark case, where one uses and 
checks the chiral theory, one uses and checks the Heavy Quark Effective Theory (HQET) when 

dealing with heavy quarks. An example is provided in fig. 3, taken from p0| . HQET predicts 

1/2 

a finite limit of fpM P , up to logarithms, for a heavy-light pseudoscalar meson P the mass of 
which M P — * 00. The two points at the extreme left in fig. 3 corresponds to a direct calculation 
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6 2 


200 + 15 + 45 


1 1 1 (3) 


APE fSl 




350 + 40 + 30 


1 14(4") 


APE (S-G) 




328 + 36 


1 1 Qf, c i N l 


UKOCD fCl 


6.2 
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UKQCD (S-C) 


6.2 


oco+lti +1U5 
Z0 "J-15 -14 


1-1415 


BLS (S) 


6.3 


235(20) ±21 


1.11 + 0.05 


Allton(S-W) 


6.0 


310 + 25 + 50 


1.09 + 0.04 


HEMCGC 


5.6 


200 ± 48 





Table 4: S refers to static quarks (infinite mass). The other values for fs have been interpolated 
to the physical B mass as shown in fig. 3. It is clear in fig. 3 and in this table that the static 
values are higher than the interpolated ones. The HEMCGC unquenched result agrees with 
those obtained in the quenched approximation. For the quenched case, (3 = 6.4 corresponds to 
a lattice spacing of ~ 0.05 fm, while (3 = 6.0 corresponds to ~ 0.1 fm. 



for infinite quark mass, according to a method proposed by Eichten |2T| which uses static quarks 
and was one of the historical starting points of HQET. The errors on these points are rather 
large. The other points correspond to finite mass, moving quarks really computed on the lattice, 
i.e. to quarks in the charm region with masses ranging up to 2.5 GeV. It appears clearly from 
the plot that all the results are consistent, and that they show that 0(1/Mp) corrections to 
the asymptotic limit are sizeable (the slope of the straight line). The 0(1/Mp) corrections are 
estimated from the curvature of the parabola. A vertical line indicates the position of the B 
meson. It may be seen that the result for is given by the interpolation between these two 
series of points. It is clear that the accuracy of the prediction for depends on the accuracy 
on fp for Mp in the Mp> region. Hence, a theoretical improvement and experimental test of 
the predictions for D decays is also useful to increase the accuracy of the predictions in the B 
sector. Results from different collaborations on fp, summarized in jjI8 |, are given in table 4. 

There exists also a specific lattice technique fit for the study of the quarkonia, and par- 
ticularly the spectra of the T's (i.e. bb bound states). This is called Non Relativistic QCD 
(NRQCD) P2| . This method is based on a discretization of the expansion of QCD in the heavy 
quark velocity: 

S NR = ^(x)DMx) + ^{x)^ + -. (16) 

where ip are two-component spinor fields (quark fields are decoupled form antiquarks). This 
method allows also to study the first hadronic excitations, while the standard method is confined 
to the ground states as exemplified in section 1270. The results seem rather good but I will skip 
a more detailed account of this field for lack of time. For recent reviews, see Hl7 . 



6 Conclusions and prospects. 

Lattice QCD is a rigorous non-perturbative calculational method based on the first principles 
of QCD. It allows to compute numerically many simple processes implying hadrons, and to 
deal with a large number of phenomenological and theoretical issues. 

The practical limitations in the computing power lead to the frequent use of the quenched 
approximation. They also yield limitations to the domain of lattice QCD: the momenta are 



quantized to a few allowed values due to the finite volume, and also bounded to values below 
or equal to ~ 1 GeV. Only the s and c quarks have masses that allow a direct study on the 
lattices. Generally speaking, only ground states and only systems with one hadron at a time 
are available. 

The light quark physics is obtained through a chiral extrapolation. Quenched approxi- 
mation becomes wrong very close to the chiral limit. The b quark physics is obtained through 
an extrapolation (or interpolation) guided by the HQET. 

The comparison with experiment is very encouraging within the present uncertainties 
which are still large: the statistical error can be as small as 5% with present statistics but the 
systematic error is seldom less than 1 5%. The unquenched calculations have still larger errors. 

The state of the art is improving quite fast, thanks to a very active theoretical work, 
and to computer improvements. Much of the latter is due to dedicated computers built by 
physicists. The Gigaflop frontier has been overcome since a few years, and the present leaders 
are between 10 and 100 Gigaflops. The Teraflop goal is within sight. The average improvement 
rate is of a factor 10 every 4 years, which means a doubling of the number of lattice points in 
every directions every 5-6 years, or a decrease of the statistical error by a factor 4 in the same 
time. The theoretical work consists in improving the lattice actions to obtain a faster approach 
to the continuum limit (a — > 0), in increasing our understanding of the finite volume effects, in 
improving the renormalization procedure, etc. In short, the theoretical improvements aim at 
lessening the systematic errors. It is of course more difficult to quantify the improvements in 
this domain than for the statistical errors. 

What accuracy will be reached in, say, ten years? I do not doubt that the statistical error 
will be as low as 1%, at least for the "best" physical quantities. Concerning the systematics, it 
is more difficult to guess. 5% ? 3% ? 1% ? One has the right to dream! Anyhow I can assure 
you of one thing: the lattice community does work hard. 
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